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^ ' Random matrix ensembles with orthogonal and unitary symmetry correspond to the 

cases of real symmetric and Hermitian random matrices respectively. We show that 
the probability density function for the corresponding spacings between consecutive 
eigenvalues can be written exactly in the Wigner surmise type form a{s)e~^^'^' for 

K* ' a simply related to a Painleve transcendent and h its anti-derivative. A formula 

p^ ' consisting of the sum of two such terms is given for the symplectic case (Hermitian 

o' 

O 

§ ; PACS numbers: 05.45.+b, 05.40.+J 
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matrices with real quaternion elements). 



It is well established that universal features of the spectrum of classically chaotic quantum 

systems are correctly described by random matrix ensembles of an appropriate symmetry [||, 

§, §. Generically there are three symmetry classes corresponding to two distinct time reversal 

symmetries plus the situation in which time reversal symmetry is absent. The level repulsion 

vj ■ — which is a characteristic of the spectra of chaotic systems and is not present in the spectra 

\^ • of integrable systems — differs for each of the symmetry classes. Thus let p{s) denote the 



probability density function for the spacing between consecutive levels, calculated after the 
energy levels are first rescaled so that the mean spacing is unity. Then p{s) oc s for a time 
reversal symmetry T such that T^ = 1, p{s) oc s^ in the absence of time reversal symmetry 
while p(s) oc s^ for a time reversal symmetry such that T^ = — 1. It is therefore convenient to 
distinguish the three cases by the label (5 and so write pp{s), where /3 = 1, 2 or 4 depending on 
the small s behaviour oipp{s). In this work succinct expressions for Pfs{s) will be given in terms 
of Painleve transcendents. 

The full probability density function Pf^{s) is by far the most studied statistic in relation 
to empirical data. For example in the early work B (this article is reprinted in [l|) on the 
energy levels of complex nuclei one sees in Figure 3 an empirical bar graph of pi{s) obtained 
from experimental data plotted on the same graph as a theoretical approximation to pi{s) 
known as the Wigner surmise. This approximation is given by the functional form Pi{s) = 
(7rs/2)e^''^*'^) . The fact that the Wigner surmise is an approximation rather than exact was 
soon realized [||, and the task of calculating the exact form oipi{s) was undertaken 0]. Actually 



the quantity calculated was £'/3(0;s), the probability there are no eigenvalues in an interval of 
length s. One then computes p/3{s) using numerical differentiation via the formula 

Pf,{s) = -^Ep{0;s). (1) 

The quantity -E'i(0; s) was tabulated by first obtaining a formula involving the eigenfunctions and 
eigenvalues of a certain integral operator. The same general approach can be taken to compute 
p^(s) for P = 2 and /? = 4 [^. The work Q revealed that the Wigner surmise, although an 
excellent approximation, can be in error by up to 2% from the exact value. 

There is a celebrated example of the empirical determination of p2{s) which is so accurate 
that the exact value to an accuracy of three to four decimal places is essential. This occurs 
in Odlyzko's numerical computation ||^, ^ of the large zeros of the Riemann zeta function on 
the critical line, in particular zero number 10^'^ and 10^ of its neighbours, which according to 
the Montgomery-Odlyzko law [^ are conjectured to have the same statistical properties as the 
eigenvalues of large dimensional random matrices. 

Nearly 20 years after the numerical computation of Ei{0; s), Jimbo et al. [^ computed the 
exact functional form of £'2(0; s) in terms of a Painleve V transcendent. Thus with a{s) defined 
as the solution of the nonlinear equation 

{sa"f + A{sa' - a){sa' -a + (a'f) = 

subject to the boundary condition a{s) ~ — s/vr — (s/vr)^ as s ^ 0, it was shown 

E2{0;s) = e^p( r^dt). (2) 



Furthermore, using a known inter-relationship between E2 and Ei [11| a formula equivalent to 



..,„.). (..(0;,)-e.p(lf(-A^)-..) ,3) 



was presented. With E2 and Ei determined, E4 can be computed from the formula [1^ 



E.(0;./2) = i(i.,(0;.) + ||i|), (4) 

Thus the exact functional form of p/3(s) can be obtained by substituting (g), (|3|) and @ as 
appropriate in (ffl) and computing the second derivative. However the resulting formulas lack 
the aesthetic appeal of (|2|)-(^), and from a practical viewpoint have the drawback of requiring 
not only the computation of crix) but also its first and second derivative. 

In this work we will show that for /? = 1 and (3 = 2 the derivative 

d „ , 

can be written in a form analogous to (^), thus allowing expressions for pp{s) to be obtained 
which have the Wigner surmise type structure ppis) = a(s)e~ '^^ A similar result will be 
obtained in the case /3 = 4. The starting point for our calculation in the case /3 = 2 is the 



formula @, but in the case /? = 1 we use a formula distinct from (y). This latter formula 
requires introducing a Painleve transcendent as (the use of the subscript B is motivated by the 
relation of this function to the Bessel kernel ||l^ ) which satisfies the nonlinear equation 

{sa'^f + a'sian - sa'^^ia's - 1) - ]{a'j,f = (5) 



subject to the boundary condition 



aB{sj r h -^. (6) 



In terms of this function we have [0] 

/ /•(W2)2 (jr(x) \ 

E^{0;s)=exp(- ^^dxj. (7) 

Our objective is to express the derivatives of (||) and (^) in terms of functional forms anal- 
ogous to the original expressions. This is achieved by using some mathematical results [15] 
relating to certain second order differential equations, of the second degree in y" , which possess 
the Painleve property. We find the functional forms involve different Painleve transcendents to 
those occuring in (^) and (j^. In the case /3 = 2, the required Painleve transcendent is specified 
by the solution of the nonlinear equation 

s\a"f + 4{sa' - a){sa' -a + (a'f) - A{a'f = (8) 

subject to the boundary condition 

For (3=1 the corresponding Painleve transcendent is specified by the solution of the nonlinear 
equation 

s\a'^f = (4(a^)2 - a'j,){s~a's - as) + -^{a'^f -^'b + I (9) 

subject to the boundary condition 

In terms of the transcendents specified by (^) and (^) our results are 

Recalling (^), (|^) and (||) we therefore have 

crfvrs) r" ait) , , , 

P2{s) = ^^exp/ -^dt 13 

s Jo t 

pi{s) = exp - / — - — dt). (14) 
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In the case /3 = 4 we start with the formula |14] 

E4{0;s) = -exp(- -J!^l dx) + - exp ( - -^^±^^ dx) 15 

where 0"^+ satisfies the same d.e. (|5|) as us, but is subject to the boundary condition 

„3/2 9/1\2/ 

.,,(.). _(i + o(<,)) + -(-)«'(i + o(,: 

Proceeding as in the derivation of ( p!2[ ) we can derive the result 



where a"_B+ satisfies the differential equation 

»'(*B+)' = (ll^W)' - Sb+Xs^W - *B+) + y(*«+) - 2*«+ + 4 C^' 



*B+w ~ 1(1 +ow) + -3^(1 + o(.: 



subject to the boundary condition 

We then have 

V,{s) = 2pi(2.) + ^ (a5+((^.)^) - 1) exp (^ - |^^"^' ^^ dt) . (18) 

The formulas (|l^), (|l^ and (15 are our main results, giving exact functional forms of a 
Wigner surmise type structure for the universal spacing probabilities pi(s), p^i^s) and ^4(5) — 
2pi(2s). These expressions are well suited to the numerical tabulations of the ppis), or the 
generation of power series expansions thereof, although that is not our concern here (accurate 
tabulations can be found in IIigII). 



Let us now turn to the derivation of (p^), (O) and (|T8|. Consider for example (|T3) (the 
derivation of ([l^) and ( [l^ ) is similar; also the derivation of ( p!3| ) will be presented as part of a 



forthcoming publication |T^). The key ingredient is the fact from |15| that the second order 
second degree equation 

x\y"f = -A{y'f{xy' - y) + A^ixy' - y) + A^y' + A^ (19) 

is solved in terms of a particular Painleve V transcendent u satisfying 

u"=(:^ + ^)M^-^ + ^^(a" + ^)+2^. (20) 

V2nn— 1/ X x^V u' x 

This is the Painleve V equation with 5 = in standard notation. As an aside we remark that 
it is known [111] that (|20) can always be solved in terms of a Painleve III transcedent. The 
relationship between (O) and (20) is via the formulas 

1 / xv! \ 2 1 — r, ^ B u — 1 'yx u + 1 ,, 

y = — ( --u) --{l-V2^)^(u-l)-^ + -^ (21) 

^ Au\u-1 / 4^ ^ ^ ' 2 u 4u-l ^ ' 



. X ( , I — u{u — 1)\2 (3 u — 1 7 

y = 



4u{u - 1) V X / 2x n 4 ^ ^ 

2 ,1 , „,2 



4 ' " ' V ' 2 



^3 = 7(/3+ ^(1 - V2^?), ^4 = ^( - /?+ ^(1 - V2^?). (23) 
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Now it is easy to see that (|19D reduces to 



X 



A. 



1 
16' 



A. 



Substituting (p^) in (p3|) gives 



V2a = 1, 



/9 



if we 


write 


(o-fi - 


1 M 


1 
'l6' 


^^ = 128 


1 

"8' 


1 

^=2' 



while use of (p3) and (p6|) in (ETI) and (22) allows us to deduce 



a 



B 



u 



1 



(TB S 

Furthermore we observe from (pl|), (p2| ) with the substitutions ( p4[ ) and (25) that 

1 



CJB + (U - 1) + 



CTB 



is also of the form ( pi|) , ( ^2|) but with 



16' 



V^=-l, /3 



(24) 
(25) 

(26) 
(27) 
(28) 
(29) 



These last three values substituted in (^) gives A2 = j^-, A^ = ^^ A^ 



^, and these values 



together with the first two identifications in ( p9|) substituted in (|T9|) give the differential equation 
(pi). The equation (O) follows from ( pTJ ) and (28), together with the facts deducible from (0) 
that as ~ as s ^ 0, while 7ro"B(s)/s^'^ ~ 1 and the boundary condition ( [lO| ) follows similarly. 

We remark that the same procedure starting with as instead of ctb does not lead to a simple 
formula analogous to (pTJ), so we cannot expect the derivative of the RHS of ([l^) to simplify; the 
results (13), (PH) and (|l^) appear to be the simplest functional forms possible. We emphasize 
that the structure of these exact functional forms for pi(s), P2{s) and Pa{s) — 2pi{2s) are of the 
Wigner surmise type a(s)e~^'*', where instead of a{s) and h{s) being simple power functions as 
in the approximation of Wigner, a{s) and b[s) are Painleve transcendents. 
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